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Abstract. We construct the iV-solitons solution in the Novikov- Veselov equation from the 
extended Moutard transformation and the Pfaffian structure. Also, the corresponding wave 
functions are obtained explicitly. As a result, the property characterizing the A^-solitons 
wave function is proved using the Pfaffian expansion. This property corresponding to the 
discrete scattering data for A^-solitons solution is obtained in [arXiv:0912.2155] from the 
9-dressing method. 
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1 Introduction 

The Novikov- Veselov equation [3, 9, 34, 41] is defined by 

U t = d 3 z U + B Z U + 3B Z (VU) + 3B Z (V*U), (1) 
B z v = d z U, d z V* = B z u. 

When z = z = x, we get the famous KdV equation (U = U = V = V) 

U t = 2U XXX + 12UU X . 
The equation (1) can be represented as the form of Manakov's triad [24] 

H t = [A, H] + BH, 
where H is the two-dimension Schrodinger operator 

h = d z d z + u 

and 

A = d 3 z + vd z + B 3 z + vB z , b = v z + v- z . 

It is equivalent to the linear representation 

H^ = 0, dti> = Aip. (2) 

We see that the Novikov- Veselov equation (1) preserves a class of the purely potential self- 
adjoint operators H. Here the pure potential means H has no external electric and magnetic 
fields. The periodic inverse spectral problem for the two-dimensional Schrodinger operator H 
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was investigated in terms of the Riemann surfaces with some group of involutions and the 
corresponding Prym ©-functions [5, 10, 22, 27, 28, 33, 37]. On the other hand, it is known 
that the Novikov-Veselov hierarchy is a special reduction of the two-component BKP hierarchy 
[23, 36, 40] (and references therein). In [23], the authors showed that the Drinfeld-Sokolov 
hierarchy of D-type is a reduction of the two-component BKP hierarchy using two different 
types of pseudo-differential operators, which is different from Shiota's point of view [37]. Also, 
in [26], it is shown that the Tzitzeica equation is a stationary symmetry of the Novikov-Veselov 
equation. Finally, it is worthwhile to notice that the Novikov-Veselov equation (1) is a special 
reduction of the Davey-Stewartson equation [20, 21]. 

Let Hip = Hlj = 0. Then via the Moutard transformation [1, 29, 30, 31] 

U(z,z) — >• U(z,z) = U(z,z) + 2ddlnu, 

ip — > 8 = — J (ipdu — ujdtp)dz — (t/jBuj — u>8ifj)dz, (3) 

one can construct a new Schrodinger operator H = d z d z + U and HQ = 0. We remark that 
the Moutard transformation (3) is utilized to construct the TV-solitons solutions of the Tzitzeica 
equation [15]. 

The extended Moutard transformation was established such that U(t,z,z) and V(t, z, z) 
defined by [13, 25] 

U(t,z,z) = U(t,z,z) + 2dd In W{?P,lo), V(t,z,z) = V{t, z, z) + 2<9dln W{ip, u), 
where the skew product (alternating bilinear form) W is defined by 

W(tp, u) = J (ipdoj - ud^)dz - (ifjdu - ujd^)dz + [ipd 3 u - ood 3 ^ + ujB 3 - ipd 3 u 

+ 2(d 2 ipdu - d^d 2 u) - 2(d 2 ipdu - Bipd 2 co) + 3V{ipdu - udip) 

-3V(ipdio -Ludip)]dt, (4) 

will also satisfy the Novikov-Veselov equation. 

In [2, 6, 7, 8], the rational solutions and line solitons of the Novikov-Veselov equation (1) are 
constructed by the 9-dressing method. To get these kinds of solutions, the scattering datum 
have to be delta-type and the reality of U also puts some extra constraints on them. In [39], 
the singular rational solutions are obtained using the extended Moutard transformation (4); 
however, the non-singular rational solutions are constructed in [4]. 

Next, we construct Pfaffian-type solutions. Given any N wave functions ^i,ip2,ip3, ■ ■ ■ ,ipN 
(or their linear combinations) of (2) for fixed potential U(z, z, t), the TV-step extended Moutard 
transformation can be obtained in the Pfaffian [1, 31] (also see [12, 35]) 



even, 



fpf^i,^,^,...,^) if JV 
\Pf(^ 2 ,</>3,---,^v) ifiVodd, 



Pf(^i, ^2,^3, • • • An) = ^€{(y)W aia2 W a3aA - ■ ■W aN _ 2UN _ 1 i) (TN , (6) 

where W ai<Jj = W{ijj a u\,il) a (j\) is defined by the skew product (4). The summations a in (5) 
and (6) run from over the permutations of {1, 2, 3, ... , N} such that u\ < 02, 03 < 04, <j 5 < erg, 
. . . and o\ < 03 < 05 < (77 < • • • , with e(<r) = 1 for the even permutations and e(<r) = — 1 for 
the odd permutations. Then the solution U and V can be expressed as [1] 

U = U + 2dd[\nP{^ u i> 2 , ip N )}, V = V + 2dd[ln P{$ 1} ifo, ^3, • • • , M], 
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and the corresponding wave function is 

^ = P(lpl,fa,4>3, • • • ,^jV,#) 
P(fa,fa,fa,...,^N) 

where $ is an arbitrary wave function different from fa, fa, fa, . . . , ifj^. 

The paper is organized as follows. In Section 2, we obtain the iV-solitons solutions using 
the extended Moutard transformation and the Pfaffian expansion. Several examples are given. 
In Section 3, the iV-solitonic wave function is derived using (7) and the Pfaffian expansion. 
Section 4 is used to prove a special property to characterize the iV-solitons wave function. 
Section 5 is devoted to the concluding remarks. 



2 iV-solitons solutions 

In this section, one uses successive iterations of the extended Moutard transformation (4) to 
construct iV-solitons solutions. 

To obtain the iV-solitons solutions, we assume that V = in (1) and recall that dd = |A. 
One considers U = -e / 0, i.e., 

dd(p = eip, tpt = ip zzz + (p Z zz, (8) 

where e is non-zero real constant. The general solution of (8) can be expressed as 

<p{z,z,t) = J e lX < lA ) v{X)dX, (9) 

where v{X) is an arbitrary distribution and T is an arbitrary path of integration such that the 
r.h.s. of (9) is well defined. 

Next, using (5) and (9), one can construct the iV-solitons solutions. Let's take v m {X) = 
<5(A — p m ) and v n (X) = a n 5{\ — q n ), where p m , a n , q n are complex numbers. Then one defines 

, _ fipm) _ J_ F(p m ) , _ ^n) _ «n 

^ m ~ V3 " V3 ' ^ ~ n V3 " y/S ' 

where 

a" ' (iA) 3 ' 

Then a direct calculation of the extended Moutard transformation (4) can yield 

W{4> m ,fa) = ia n q -^^e F ^ +F ^\ W{4> m A n ) = i^LZ^e^+'K 

Qn ~r Pm Pn T Pm 

W(^ m ,fa) = ia m a n q -^^e F ^ +F ^\ (10) 

Qn ~r Qm 

The iV-solitons solutions are defined by 

U(z, z, t) = — e + 2<9(9 In tat (z, z, t), V(z, z, t) = 2ddhiT]^{z, z, t), 
and then 

U(z, z, t) — > — e as zz — > oo, 



F(X) = (i\)z + (iXft + ~z + T^vT*- 
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where t is fixed. The r-functions are defined as follows. For simplicity, let's denote 

W(pm,q n ) = W((j) m ,tpn), W(p m ,Pn) = W((j) m , <f) n ), W(q m ,q n ) = W(ip m ,ip n ), (11) 

and notice that F(—X) = —F{\). The tn is defined as 

t n (z, z, t) = Pf(-pi,qi, -p 2 , <?2, -P3, 93, • • • , Pn i 9iv), (12) 

where 

(-Pm,-Pn) = W(-p m ,-pn), ("Pm, 9n) = W(-p m ,q n ) + 6 mn , 

(q m , q n ) = W(q m , q n ). (13) 
To get the expansion of (12), we use the following useful formula [14, 38] 



Pf(A + B) = £ E (-l) H " r Pf(A a )Pf(B aC ), 



(14) 



where A and B are m x m matrices and s = [m/2] is the integer part of m/2; moreover, we 
denote by a c the complementary set of a in the subset {1, 2, 3, . . . , m} which is arranged in 



increasing order, 


and a = 


at + «2 + • 


has 
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-1 



{-pi,qi) {-pi,-P2) (-01,92) 

(91,-02) (91,92) 

(~P2,92) 

(9JV,9l) {(IN,-P2) (9iV,92) 



(-0i,9iv) 
(9i,9Jv) 
(-P2,9at) 





where A^r and B^r are 2iV x 2iV matrices. Hence by (14) one can have the expansion of (12), 
i.e., 



N 



r N = l + Y,fi+Y,\ E fiJh---fi m n pv»i- 

£=1 m=2 \l<£ 1 <£ 2 <---<^m<A r l<j<fc<m 



(15) 



where 



9£ -0£ 



(p<,- -pikKqtj - gijfaj + 94) (ggj +pe k ) 

{Pij + Pl k ){<li j + Qijlptj ~ qe k ){<ttj ~ Pe k ) 



Here we have utilized the formula that if C is a 2N x 2N matrix with (i, j)-th entry , then 



one has the Schur identity [32, 36] 



pf(c) = n 



l<j<j<27V v J 



«i + a, 



(16) 
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Next, we illustrate the formula (15) (or 12) with several examples. 
(1) One-soliton solution: 



Ai(z,z,t) 



(-Pl,?i) 
(<?i,-Pi) 



B, 



1 

-1 



Then 



n = Pf(Ai +Bi) = 1 + iai ^±^e F ^- F ^\ 

91 - pi 



(2) Two-solitons solution: 



A 2 (z,z,t) 









(9i, -Pi) 

(-P2,-Pi) (~P2,9i) 

. (?2,-Pl) (?2,9l) 



(-Pi,9l) (-PU-P2) {-Puqz}' 



Bo 



10 
-10 
1 
0-10 



{qi, -P2) 



(?2, -P2) 



(91,92) 
(-P2,92) 





Then 



or 



r 2 = Pf (A 2 + B 2 ) = 1 + iai ^±^e F ^- F ^ + ia2 ^±® e F (»)- F ^ 

9i - Pi 92 - P2 

I ia^a / 1 +qiP2 + q2 P2 ~ Pi 92 - 9i Pi + 92 P2 + 91 gFfaO-FfrO+FfaQ-Ffra) 

91 - Pi 92 - P2 P2 + Pi 92 + gi 92 - PI P2 - 91 



T2 = l + /l + /2+Pl2/l/2, 



(17) 



where 



12 



91 - Pi 

Pi - P2 91 - 92 Pi + 92 91 + P2 
Pi + P2 91 + 92 Pi - 92 91 - P2 ' 



f 2 = ia 2 P2+q2 e F ^- F ^\ 



92 -P2 



The t 2 soliton (17) is also found in [2] using the 3-dressing method. 
(3) Three-solitons solutin: 



A 3 (z,z,t) 
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!, -P2) 







(~P3,-Pl) (~P3,9l) (~P3,-P2) (~P3,92) 
, (93, "Pi) (93,9l) (93,-P2) (93,92) 



(92, -P3) 


(93, -P3) 
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(~P3,93) 
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Then 



r 3 = Pf (A3 + B 3 ) 

= 1 + h + h + / 3 + P12/1/2 + P13/1/3 + P23/2/3 + P12P13P23/1/2/3, 



where 



h = iai Pl + qi e F M- F ^\ h = ia 2 P2 + Q2 e F ^- F ^\ 
91 - pi 92 - P2 



12 



Pi - P2 qx - q2 pi + q2q\+ P2 



"13 



93 -P3 p\+ P2 qi + q2Pi- q2 qi - P2 

pi + q3 qi + P3 pi -pzqi-q-i m P2 + 93 92 +P3P2- P3 92 - 93 



pi - 93 qi - P3 pi +P3qi + 93 
The wave functions 



23 



P2 ~ q3 q2 - P3 P2 + P3 92 + 93 



(18) 



In this section, one uses (7) to construct the corresponding wave function of the r function (15). 

From (7), one knows that the corresponding wave function of the A^-solitons (15) can be 
written as 



p( ¥(-Pi) fill) f{~P2) p(q2) 
\ V3 ' V3 > V3 ' V3 : 



TN 



Using the notations in (11), (12) and (13), we can express cpw as 
P(~P1, 91, ~P2, 92, • • • , ~PN, 9Af, A) 



But we notice that 

(-p m , A)" = W{-p m , A), (9m, A)" = W{q m , A), 



(19) 



(20) 



where (•)" means there is no 5 mn here when compared with (13). Now, let's compute P{— pi, 9i, 
— P2-, 92, • • • , —Pn, Qn, A) using (14). In this case, 



P(~Pl,qi, -P2,92, 



~Pn, 9w, A) = Pf (Mat + Q 



N, 



where 



M N 



A N (z,z,t) 



(A,-pi) (A,9i) 
y(-pi) y(gi) 



and 



A? 








y/3 















(A,9a0 

<p(<in) 

^3 



{.-pi A) 

(91, A) 

(9AT,A) 



V3 



f(-pi) ' 

V3 

y(gi) 



V3 

y(A) 

^3 
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Using (16), a simple calculation yields 
Pf(M 7V + Qn) 



N 



N 



he, he,, • • • he 



n 



m=2 \i£e 1 <e 2 <£3<-<e m SN 



l^j<k^m 



0Xjv(A), 



(21) 



where 



¥>(A) 



h e (X) = iat 



Pe + qePe + Xq e - X 



\/3 ' " "q e -piPi- Xq e + X' 

N N 



,F(q e )-F(p e ) 



X N (X) = l + ^hz(X) + Y^ 



E 



he, he„ • • • hp 



t=i 



m=2 \l^ 1 <£ 2 <£ 3 <...<£ m ^Ar 



Pi + Xqi-X 
Pi - X qi + A ' 



and IV 4 is defined in (15). 

We give several examples here. 
(1) The one-soliton wave function: 



Mi 



(-Pi,A) 
(ft, A) 

(A,-pi) (A, ft ) 

y(-pi) y(A) 
V3 ^ 



^i(z, 



V3 

y(gi) 
y(A) 

^3 




Qi 



Bi 













Then 



P(- Pl , qi ,X) Pf(Mi + Qi) 



1 + . Qi Pi+9iPi + Agi-A cF((;i) _ Ffa) 
qi - Pi Pi - A qi + A 



We remark that 

¥>1 



l + /i 



n . ■ Pi + A qi - A 

1 + mi — /i 

pi - X qi + X 



1 + /: 



1 + iai 



2pi 



Pi - A 



2?i 



qi + A 



l /i 



(l + /i) + 2*oi( 



pi 

Pi— A qi+Av 



_23_V F (9i)-^(Pi) 



1 — 2iai 



l + /l 



Pi 9i 



A - pi A + qi 



,F{ qi )-F{ P1 ) 



Tl 



(22) 



This is the one-soliton wave function in [2, p. 9]. 
(2) The two-soliton wave function: 



Mo 



A 2 (z,z,t) 



(A, — pi) (A,gi) (A, — p 2 ) (A,g 2 ) 
y(-pi) y(-pa) g(gg) 

%/3 v 7 !? 



(-J*,A) ^ 
(9i, A) 



y(gi) 

<P(-P2) 



(-ft,, A) 
(92, A) 





y(A) 
' V3 



y(gg) 
y(A) 
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Q2 



B 2 



0' 













Then using (21), one has 

P{—Px,1ii ~P2, <72, A) Pf(M 2 + Q 2 ) 



<T2 



T2 



T-2 



72 



[l + h 1 (X) + h 2 (X)+F 12 h 1 (X)h 2 {X)] 



1 + iai Pl + gi Pi + A gi ~ A c F( gl )-F( Pl ) + , a , P2 + 92 + A g 2 - A cF(g2) _ F(p2) 



?i - Pi Pi - A qi + A 



' 92 - P2 P2 ~ A q 2 + A 



T2 

I ia,ia 2 Pl + 91 yi + A Ql - X P2 + Q2 P2 + A Q2 - X c ^)+^)-F( Pl )-F( P2 ) 



' qi - Pi Pi - A qi + A q 2 - p 2 p 2 - A q 2 + A 
(3) Three-solitons wave function: 



(-Pi,A) 
(ft, A) 



A 3 (z,z,t) 



5 a 



(A, — Pi) (A,gi) (A, -p 2 ) (X,q 2 ) (A,-p 3 ) (A,</ 3 ) 
y(-pi) y(-P2) yfe) y(-P3) yfe) 

V3 \/3 a/3 \/3 \/3 %/3 

0' 











00000000 
00000000 



From (21), we get 

P{~Pi,Qi, ~P2, q 2 , ~P3, <?3, A) Pf(M 3 + Q 3 ) 



y(-pi) ' 
y(gi) 

(-P2,A) ^fl 
(ft, A) 

(-P3,A) ^f) 

(ft, A) « 

n i( A ) 

_MA1 n 
Va u 



<5s 



r 3 



T~3 T 3 

[1 + /n(A) + h 2 (\) + /i 3 (A) + Pi 2 /»i(A)/i 2 (A) + Pi3^i(A)/t 3 (A) + P 23 /i 2 (A)/ i3 (A) 



where 



*1 
^3 



+ P 12 P 13 P 23 / ll (A)/i 2 (A)/i 3 (A)] 



,^ pi + Agi - A Pl + gi c j- fa )- Ffpi ) 
pi - A qi + A - pi 

ia _ P3 + A g 3 - A p 3 + g3 e F( g3 )-F( P 3) 

P3 - A g 3 + A g 3 - p 3 



■ Q , P2 + A g 2 - A p 2 + g 2 cF(q2) _ F(P2 ) ^ 
p 2 - A g 2 + A q 2 - p 2 



and 



"12, ri 3 



and P 23 are defined in (18). 
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4 A property of iV-solitons wave function 

In this section, we will express the wave function (19) as another form to generalize the equa- 
tion (22) to iV-solitons case. 

Firstly, according to the Pfaffian expansion in [11], it is not difficult to see that 

Pf(6i, b 2 , h, h, ■ ■ ■ , b 2n -i, fon, hn+i) 

2n+l 

= E (-iy +m (bj,b m )Pi(b 1 ,b 2 , . . . , bj, . . . , b m , . . . , b 2n , b 2n+1 ), (23) 

m=l 

for j = 1,2, ...,2n + l, 

where bj and b m mean these two terms are omitted. 

Secondly, noticing (10) and letting A = — p a or A = q a , a = 1, 2, . . . , n, we have 

Pf(-Pl, Ql, ~P2, Q2, ■ ■ • , ~Pa, Qa, ~Pa+l, Qa+1, ~PN, QN , ~Pa) 

= Pf(-Pl, qi, ~P2, <72, • • • , ~Pa, -Pa+l,q a +l, ~PN, <?iv) = 4>{~Pa)XN (~Pa) , 

Pf(-pi, qi, -p 2 , q2, ■ ■ ■ , q a -i, -Pa,q a , -p a +i,qa+i, -pn, qN, q a ) 

= Pf(-pi,gi, -P2,q 2 , ■ ■ .,q a -i,q a , -p a +i,q a +i, ~PN,qN) = a a 4>{q a )xN{q a )- (24) 
They can be seen as follows. By (23), one has 

Pf(-pi, qi, -P2, q2, ■ ■ ■ , ~Pa, q a , ~Pa+i, q a +i, ~Pn, qN, ~Pa) 

= ~{qa, -pi)Pf(<?l, ~P2, q2, ■ ■ ■ , -Pa-l,qa-l, ~Pa, ~Pa+l, ~PN, qN, ~Pa) 

+ {q a ,qi)¥i{-pi, ~p 2 ,q2, ~Pa-l,q a -l, ~Pa, ~Pa+l, ~PN,qN, -Pa) 

- {qa, -Pa)¥t{-Pl,qi, ~P2,q 2 , -p a -l,q a -l, -Pa+l, ~PN,qN, ~Pa) H 

+ {q a , -p ct ) tt Pf(-pi,gi, ~P2,q2, -p a -l,q a -l, ~Pa, ~Pa+l, ~PN^n) 

= [-{Qa,-Pa) + {q a , -p a f]P£{-Pl,qi,-P2, q2, ■ ■ ■ ,-Pa-l,q a -l,-Pa, ~Pa+l, ■ ■ ■ ,~PN, <?7v) 
= Pf(-pi,<7i, ~P2,q2, -Pa-l,q a -l, ~Pa, -p a +l, ••• , ~PN, 9iv), 

where {q a , — p a )^ is defined in (20) and we know that 
Pi{. . . , -p a , . . . , -p a ) = 0. 

The second equation of (24) can be proved similarly. 
Finally, from (23), one yields 

Pf(-Pl, qi, -P2, q2, ■ • • , Pn , qN, A) = (A, -pi)Pf {qi, ~P2,Q2, ~Pn, qN) 

- (A,<?i)Pf(-pi,-p2,<?2, • • • , -pn,qn) + (A,-p2)Pf(<?i,-f>i,<?2, • • -,~PN,qN) 

- (A, q 2 )Pi{qi, -pi, ~P2, ~Pn, qN) H 

+ (A, -p N )Pi{-pi, -p 2 , q 2 , ■ ■ ■ , ~PN-i,qN-i,qN) 

- (A, q N )Pi{-pi, -p2, <72, • • • , -PN-uQN-l, ~Pn) + ^(A)^. 

Also, the wave function (21) can be written as 

ip N {\) = 4>{\)xn{\), 
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where Xiv(A) = ■ Therefore, using (24) and letting AF n = F(q n ) — F(p n ), we get 



/ x \ -i Pi + A . AF , , qi — A . A Fi 

Xat(A) = 1 - zaie 'XivWlJ r^aie ^(-pi) 

9l + A 



A-pi 



za 2 e 2 xtv(92) 



A -P2 

P^v + A . AF 



92 + A 



ra 2 e 2 Xw(-P2) 



X-p N 
2p 



ia N e n Xn(qn) 



Qn — A . AF , . 

— wjve n Xn{-Pn) 

Qn + A 

2gi 



A-pi 



1 ■ AFi / \ z 9l ■ AFi / \ 
— zaie ^tvWi) — r^aie 1 xn{-Pi) 



2p 2 

X-P2 

2p N 



ia 2 e AF2 XN(q2) 



ia N e AFN XN(qN) 



qi + X 

2 92 . AF 2 ( \ 

92 + A 
2q N 



ia N e AFN XN(~PN) 



X-PN QN + A 

+ [-iaie AFl X7v(<?i) + iaie AFl XN(-Pi) ~ ia 2 e AF2 XN{q2) 

+ ia 2 e AF2 XN(-P2) iaNe AFN XN(qN) + ia N e AFN 'xn(-Pn)] ■ 



(25) 



Since Pf (— pi, qi, — p 2 , q 2 , ■ ■ ■ , — Pn-> Qn-> 0) = ^/§> we nave Xiv(0) = T JV- Then the last term in 
[• • •] of (25) is zero (or lim Xiv(A) = 1 )• Hence one has 

A— >oo 



Xiv(A) = 1 



- t iaie 1 XN{qi) ~^iaie 1 xn{-Pi) 

X-pi qi + X 

2 P2 . af 2 i \ 2q 2 . AF 

ia 2 e 2 XN{q2) —ia 2 e 2 Xn{-P2) - 

X - p 2 q 2 + A 

2PN . AFn i \ %QN 



X-PN 



ia N e N XN(qN) 



ia N e AFN XN(~PN)- 



qN + A 



This formula is also obtained by the d-bar dressing method when the d-bar data is the degenerate 
delta kernel [2, p. 6]. 

When n=l, we have (22). For n=2, from (21), one knows that 



X2(-Pi) 
X2(<?i) = 

X2{~P2) 
X2fe) = 



p 2 +q 2 p 2 -pi q 2 +pi A F 2 
1 + ia 2 e 2 



92 - P2 P2 + Pi 92 - Pi 

. . . P2+q2P2+qiq2~qi af 2 
1 + ia 2 e 2 

92 -P2P2- 91 92 + 91 
-. . • Pi + 91 Pi ~ P2 91 + P2 AF 1 

1 + ia\ e 1 

9i - Pi Pi + P2 9i - P2 

, . ■ Pi + 91 Pi + 92 91 - 92 A F\ 
l + iai e 1 

9i - Pi Pi ~ 92 9i + 92 



A"2, 
A-2, 
A-2, 

/r 2 - 



Then 



/w i • AFi i \ 2 9i . AF , > 

X2(A) = 1 - - ia x e 1 X2{qi) -^ia\e 1 X2{-Pi) 

X-pi 9i + A 



2 P2 



X -p 2 



ia 2 e AF2 X2{q2) 



2 92 . af 2 i \ 

— ria 2 e 2 X2{-P2) 

92 + A 



We remark that this formula also appears in [2, p. 10], the parameters being different. For n=3, 
by (21), we obtain 



X3(-Pl) 



p 2 + 92 P2 ~ Pi 92 + Pi AF 2 , • P3 + 93 P3 ~ Pi 93 + Pi AF 3 

1 + ia 2 e 2 +ia% e 3 

92 - F2 P2 + Pi 92 - Pi 93 - P3 P3 + Pi 93 ~ Pi 
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P2 +q2P2~ Pi g2 +P1P3+ g3 P3 ~ Pl g3 + Pl ro AF2+AF3 

. . r 23e 



q2 -P2P2+ Pl q2 - Pl <?3 - P3 P3 + Pi g3 



xs(qi 



P2 + q2P2 + qiq2- qi af 2 , • P3 + g3 P2 

l + ia 2 ; e 2 +ta 3 



' <?2 - P2 P2 - qi q2 + qi 



g3 - P3 P2 



Pl 

gl g3 ~ gl c AF 3 



/r 3 , 



. . . P2 + g2 P2 + gl g2 - gl P3 + g 3 P2 + gl g3 ~ gl TO AF 2 +AF 3 

+ lCt2ia3 p ""° 

g2 -P2P2- gl g2 + gl g3 - P3 P2 - gl g3 + gl 



gl g3 + gl 

P 23 e z 



X3(-P2) = 



1 + icq 



Pl + gl Pl - P2 gl + P2 AFi , • P3 + g3 P3 - P2 g3 + P2 AF 3 

e 1 + 203 e 3 



gl - Pl Pl + P2 gl - P2 



P3 - g3 P3 + P2 g3 - P2 



+ iaiia.3 



Pl +glPl -P2 gl +P2P3 +g3P3 -P2 g3 +P2 m AF1+AF3 

ri 3 e 



X3(g2 



X3(-P3) 



gl - Pl Pl + P2 gl -P2P3- g3 P3 + P2 g3 - P2 
Pl + gl Pl + g2 gl - g2 AF 1 , ■ P3 + g3 P3 + g2 g3 - g2 AF 3 

l+icbi e + ««3 e 3 

gl - Pl Pl - g2 gl + g2 g3 - P3 P3 - g2 g3 + g2 

Pl + gl Pl + g2 gl ~ g2 P3 + g3 P3 + g2 g3 



+ ia\ia 3 - 



92p ige AFi+AF 3 



/T3, 



1 + ia\ 



■ Pl Pl - g2 gl + g2 g3 - P3 P3 - g2 g3 + g2 
Pl + gl Pl - P3 gl + P3 AFi , • P2 + g2 P2 - P3 g2 + P3 AFo 

e 1 + id2 e 2 



+ ia\icL2 



Pi ~ gi Pi + P3 gi - P3 
Pl + gl Pl ~ P3 gl + P3 P2 



X3(g3) = 



Pl - gl Pl + P3 gl - P3 g2 ~ 
, . - Pl + gl Pl + g3 gl - g3 AFi , 

1 + ia\ e 1 + ia 2 

gi - Pi Pi - g3 gi + g3 g2 



g2 - P2 P2 + P3 g2 - P3 
+ g2 P2 ~P3 g2 +P3 TO AF1+AF2 

ip 12 e 



P2P2 +P3 g2 -P3 

P2 +g2 P2 +g3 g2 ~ g3 AF 2 
P2 P2 - g3 g2 + g3 



/T3, 



+ ia\ia2 



Pl + gl Pl + g3 gl ~ g3 P2 + g2 P2 + g3 g2 ~ g3 p g AFi+AF 2 
gl - Pl Pl - g3 gl + g3 g2 - P2 P2 - g3 g2 + g3 



/T3, 



where P12, P13 and P23 are defined in (18). Then 
2 Vl ._ _AFi.. \ ^g 

2p 3 



X3(A) 



1 



1 ■ AFi t \ z gl • AFi / \ 

— ^ale 'xslgi) —riaie ^(-pi) 

A - Pi gi + A 



2p 2 
A -P2 



ia,2e 



AF 2 



X3(g2 



2g 2 
g2 + A 



ia>2e 



AF 2 



X3(-P2) 



ia 3 e Ai?3 X3(g3) 



A -P3 



2g 3 
g3 + A 



ia 3 e AF3 X3(-P3)- 



5 Concluding remarks 

In this paper, we have used the extended Moutard transformation to construct the iV-solitons 
solutions. The basic idea comes from the successive iterations of solitons solutions, as remains to 
be the simple method to obtain the iV-solitons solutions. Also, the corresponding wave functions 
are constructed by the Pfaffian expansion of the sum of two anti-symmetric matrices (14) when 
compared with the 9-dressing method [2]. 

To obtain real iV-solitons solutions of the Novikov-Veselov equation (1), one has to put 
extra relations between — pi and qi [2] . It could be interesting to investigate these real solutions 
for the Schrodinger operator (self-adjoint). On the other hand, the resonance of iV-solitons 
solutions of DKP or KP theory has been studied in [16, 17, 18, 19]. And then the resonance 
of iV-solitons solutions of Pfaffian type (15) deserves to be investigated. These issues will be 
published elsewhere. 
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